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Abstract 

We study the density perturbation by a curvaton with a double well potential and estimate 
the nonlinear parameters for non-Gaussianity and the amplitude of gravitational wave background 
generated during inflation. The predicted nonlinear parameters strongly depend on the size of 
a curvaton self-coupling constant as well as the reheating temperature after inflation for a given 
initial amplitude of the curvaton. The difference from usual massive self-interacting curvaton is 
also emphasized. 
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I. INTRODUCTION 



Cosmic inflation solves various problems in the standard Big Bang cosmology such as 
flatness, the horizon, and the monopole problems Simultaneously, the quantum fluctua- 
tion of a light scalar field, e.g., inflaton field 0, generated during inflation is stretched by the 
rapid cosmic expansion and provides the seed of large scale structure in our Universe [2J. 

The density perturbation generated in a single field inflation model is scale-invariant and 
almost Gaussian with the corresponding nonlinearity parameter /nl much less than unity . 
This is consistent with the current limit on the local type non-linearity parameter /nl from 
the Wilkinson Microwave Anisotropy Probe (WMAP) seven-year data, —10 < /nl < 74 



at the 95% confidence level 



4j. It is expected that the observational sensitivity is gong to 



improve significantly by the Planck data j^j and using large scale structure data within the 
near future. The non-Gaussianity could be an important observable to discriminate between 
various mechanisms of density perturbation generation. 

On the other hand, beyond the simple canonical single field slow-roll inflation, the large 
non-Gaussianity with different shape are generally predicted 6|, |7J. There are many models 
for the generation of the observed density perturbation and a large non-Gaussianity. This 

15 1, preheating or deep 



can happen during inflation 
in the radiation dominated era 
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12| 



at the end of inflation 
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The last case includes the "curvaton" scenario |18H22j]. where the scalar field is too light 
to make effects around the inflationary epoch but it might play an important role much 
later in the early Universe. As the Universe expands, the cosmic expansion rate H becomes 
comparable to the mass of the curvaton and the curvaton field starts to oscillate in the 
radiation dominated era. After that when the expansion becomes less than the decay rate of 
the curvaton, it decays to light fields and the isocurvature perturbation of the curvaton field 
becomes adiabatic or mixed with that from the inflaton field. If the curvaton energy density 



is subdominant at its decay time, the large non-Gaussianity is generated in general 



Another important measure is the gravitational wave background produced during in- 



flation 



24J parametrized by the tensor-to-scalar ratio r T , because it could directly indicate 



the energy scale of inflation, tt is related with one of the inflaton's slow roll parameter e 
as tt = 16e, while the density perturbation V( oc e~ l in a single field inflation model. On 
the other hand, in the curvaton scenario, the density (scalar) perturbation comes from the 
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curvaton. Nevertheless the (non- ) observation of Tt gives strong constraint on the 



eters of the curvaton scenario 



param- 



25] . The present bound is tt < 0.36 (95% CL) [4J, which 
is expected to be tightened as ~ 10 _1 from Planck satellite Q| and ~ 10 -3 from 
DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) (26] and/or the Big 
Bang Observatory(BBO) Q. 

Curvaton scenarios have been often modeled by a scalar field a with a quadratic poten- 
tial V = \m 1 a a 1 . In that case the non-Gaussianity has been studied with the sudden d ecay 



approximation |28H30j . which shows good agreement with the full numerical approach [31]. 
Beyond the simplest model of the curvaton, there are various possibilities; the inflaton per 
turbation may not be negligible 32( , the curvaton can have different types of potential 



33 



and there could be multiple curvaton fields 



due to non-quadratic terms can be seen 
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Significant effects on non-Gaussianity 



So far, it has been assumed that the mass squared at the origin of field is positive as above. 
However, there is no reason that the true minimum is located at the origin. Scalar fields 
have been often introduced for spontaneous symmetry breaking in particle physics models. 
A moduli field, which is a promising candidate of curvaton 21], or the Peccei-Quinn field, 
to solve the strong CP problem, usually has the large vacuum expectation value (VEV). 

In this paper, we examine a curvaton model with a double well potential where it develops 
nonvanishing VEV. Since a curvaton needs to develop large expectation value during infla- 
tion, throughout this work, we assume that the potential is very flat and the self-coupling 
constant is small enough. Such a tiny self-coupling scalar field model has been studied in 

nn 

the axion model to solve the isocurvature perturbation and domain wall problems [461 . 1471 ] . 
Due to the flatness the curvaton has a large initial amplitude, cr* > v, and almost stays 
there during inflation with negligible movement. After inflation the curvaton field starts to 
roll down into the minimum, oscillates and decay into radiation. After the curvaton decay, 
its isocurvature perturbation is transferred to the adiabatic curvature perturbation in the 
radiation dominated plasma. 

We take account of the inflaton perturbation as well as that of the curvaton. Therefore 



the curvature perturbation shows generalised mixed inflaton-curvaton type [32J. We con- 
sider the initial amplitude of the curvaton field is arbitrary but larger than the symmetry 
breaking scale. In the opposite case where the field starts rolling around hill of the potential, 
interesting results have been shown 
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The paper is organized as follows. After describing the model with a double well potential 
in section [III we consider the cosmological evolution of the scalar field in section IIHI In 
section |IV] we estimate the non-Gaussianity and gravitational wave background of this model. 
We summarize our results in section [V] 



II. A CURVATON MODEL WITH A DOUBLE WELL POTENTIAL 

We consider a real scalar curvaton model with a double well potential. The Lagrangian 
density of the field is given by 

£ = - 1 -(daf-V(a) 1 (1) 

V{a) = \(o 2 -v 2 )\ (2) 

with A and v being respectively the self coupling constant and the VEV. When the field has 
a nontrivial expectation value in the potential Eq. (|2J), the effective mass of it is expressed 
as 

V aa = X(3a 2 -v 2 ) (3) 
and the mass at the true vacuum is given by 

ml = V aa \ a=±v = 2Xv 2 . (4) 

The decay rate of a depends on its interaction with light particles. If a couples with a light 
fermion ip through a Yukawa interaction as £j nt = yipaip, the decay rate is roughly given by 

v 2 

T CT ~ f-m CT . (5) 

If a does not directly couple with light particles, the decay rate would be expressed as 

T a = C^f = C{2Xf' 2 v. (6) 

with C being a numerical coefficient of including coupling constants and phase volume. This 
kind of decay rate formula is realized, for instance, for the radial direction of Peccei-Quinn 



field in the hadronic (KSVZ) axion model 49( 1 . 



1 The corresponding curvaton scenario with Peccei-Quinn field in the extension of MSSM has been studied 
in Ref. 0. 
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III. COSMOLOGICAL EVOLUTION OF a 



By definition, the curvaton field is subdominant during inflation. Its contribution can be 
important after inflation in the deep radiation dominated era. For this purpose, we consider 
the case that A is very tiny. In such a case, the potential is very flat, like chaotic inflation 
with quartic potential, and hence the fields can develop a large expectation value during 
inflation. 

The equation of motion for the homogeneous part of a is given by 

a + 3H& + X(a 2 - v 2 )a = 0. (7) 

Before the inflaton decay, the energy density of the Universe is dominated by the inflaton 
whose equation of motion is given by 

, dV 

+ 3#0+ — = 0, (8) 
and then the Friedmann equation is 

3M 2 H 2 = p^ + p CT , (9) 

with ^> p a . Here Mp ~ 2.4 x 10 18 GeV denotes the reduced Planck mass. After the 
inflaton decay, the Universe is dominated by the radiation generated from the inflaton decay, 
Then the field equations, instead, are 

p r + AHp r = 0, 

(10) 

3M 2 H 2 = p r + p a . 

The curvaton field with a large expectation value during inflation almost stays there until 
the Hubble parameter H becomes comparable with the effective mass, i.e. 

#os^K.| os = A(3a 2 s -*, 2 ). (11) 

From that time, the curvaton starts to oscillate with the initial amplitude <r os and the energy 
density 

P*L = \{°l-v 2 ) 2 . (12) 

When the initial amplitude of the curvaton is much larger than the location of the minimum 
as cr* ^> v, the evolution of the curvaton is dominated by the quartic potential and the hill 



at the origin can be ignored during the oscillation. In this case, the energy density of the 
curvaton decreased as a(t) -4 after the oscillation starts. 

After the oscillation amplitude decreases enough so that the field cannot go across the 
potential hill around the origin, the field can settle down at one of the two degenerate and 
distinct vacua. Once the curvaton field find one of the minimum, the oscillation amplitude 
becomes less than of the order of v, and the energy density becomes smaller than 

Pa\ v =\tf-f?, (13) 

with a v being the amplitude of the order of v that denotes the transition from quartic 
oscillation to quadratic one occurs. Note that a v is independent of cr*. The ratio of the 
energy densities given by Eqs. (112]) and ( TT3|) is scaled by (a v /a os )~' i , because the quartic 
term initially dominates. The precise value of a v does not affect the differentiation by cr* 
but it may affect the energy density of the curvaton when it decay, since the exact transition 
epoch between quadratic and quartic affect the evolution of the curvaton energy density. 

After the curvaton find a minimum at (a) = v or (a) = —v, its energy density deceases 
as a pressureless matter oc a(t)~ 3 since the quadratic potential dominates. Which of the 
VEV would be realized depends on the initial field expectation value. The dependence is 
shown as the function defined by 0(<t O s) = (&)/v in figure CD 

At the late time t ^> t v in the deep oscillation period dominated by the quadratic poten- 
tial, the evolution of a can be well expressed as 

a{t) ~ vQ(a os ) + ^^3/4 sinmj. (14) 

The amplitude of the oscillation o"2os can be estimated by using the simple scaling law 
between t os and t v 36[ for radiation dominated (high Tr) and oscillating inflaton dominated 



(low Tr) at H os respectively by 



3/4 

= ) , for high T R , 
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FIG. 1: The position of VEV with varying initial amplitudes of the curvaton. The vertical axis is 
the sign of VEV, {cr)/v and the horizontal axis is cr*/v. 

Here, the time of reheating can be approximated when the Hubble parameter is similar as 
the decay rate of inflation I\ and the reheating temperature T R is estimated by 



J* 2 — — n T 4 

R ~ 3M130 9 * R - 



(17) 



We consider or T R as a free parameter. Equation (fT5|) has been sometime noted as a [30] 



or cr os [3b|J in literature. Figure [2] shows the good agreement between analytic approximated 
solutions Eq. ( IT4|) and full numerical solutions. Finally, the a field decays into radiation, 
when the Hubble parameter H becomes comparable with its decay rate H ~ T a . 

When v < a* < a v , the curvaton field starts to oscillate initially in the potential domi- 



nated by quadratic term when if 2 



m 



<l „. Therefore we find that 



cr 2os ~ (tr» - ve)(m a /2H os ) 3/i . 



(18) 



Here we summarize the conditions for the a field to be a viable candidate for curvaton. 
The curvaton is almost massless and its field value is frozen during inflation. This is expressed 
by V aa H^ nf and rewritten as 



(I) 3A(a 2 - v 2 ) « 



(19) 
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FIG. 2: The comparison between analytic solutions (purple line) and numerical solutions (blue 
line) for various initial values a os = Av (left), 10f (center), and 50t> (right). The vertical axis is - 
and the horizontal axis is m a t. These figures show that analytic solutions well describe asymptotic 
evolution of the field. 

The curvaton energy density is subdominant compared with that of the inflaton during 
inflation, which is expressed as 

^-v 2 ) 2 «3M 2 H 2 n , (20) 

This condition is automatically satisfied from Eq. ( II 9p when a* < Mp. Whether a dominates 
the energy density of the Universe at the moment of a decay depends on the reheating 
temperature after inflation Tr determined by the inflaton decay rate 1^. We will pursue 
the details for this in the following subsections. During preheating, the symmetry might 
be restored and the topological defects could be formed due to the parametric resonance 



and the large fluctuation of (5a 2 ) ^> v 2 51[- This problem can be avoided if the dynamics 
is pure classical and the initial fluctuations 5a /a is less than the change of the amplitude 
of a per one oscillation when the curvaton field settles down to one of two minima of the 



potential 47j. This condition is easily written down as 



5a # inf /27T AA H c V\v 2 /a os , . 

— < — ; 7f , \ 2i -> 

a a os A uj y/Xv 

where we have used the fact that the curvaton oscillation is dominated by quartic term 

between the end of inflation and the critical point and thus c^ s /a^ = H 2 S /H 2 . Therefore 
there is no domain wall problem for a large VEV satisfying 

(II) ^ < v. (22) 
In this paper, we consider this large symmetry breaking scale. 
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A. A high reheating temperature case 



First, we consider the case that the reheating after inflation is completed when a starts 
to oscillate, which means that the Hubble parameter at reheating Hr is larger than that at 
the beginning of the oscillation. This condition of high reheating temperature corresponds 
to H R > H os or 



Tn > 



7T 



-1 



—g, ) 3M^A(3a 



,2\2 



1/4 



The energy density of radiation produced by the inflaton decay at H os is 

Pr 



Ls = 3M^ o 2 s . 



(23) 



(24) 



The energy density of the curvaton and radiation from the inflaton decay at H ~ IV are 
given by 

3 



A 



«r CT 



and 



(25) 



(26) 



The u to radiation ratio is evaluated as 



Pa 


Xv 4 


/ a v 


Pr 




W CT 




. v^aU 2 ( 


'3Ai; 2 




36M| \ 





(27) 



for the radiation dominated Universe. 
For the case of v < a os < a v , we obtain 

X l/ Ha 2 2 )2 



1271^2(3^-^)3/4^/2 

for the radiation dominated Universe. 



(28) 



B. A low reheating temperature case 

Next, we consider the case that the inflaton still oscillates (we assume the quadratic 
oscillation) around the minimum and the reheating is not completed yet when a starts to 
oscillate, Hr < H os , which is the opposite condition of Eq. (123]) . 
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The energy density of the inflaton at H os is 



3M 2 P H 2 os . 



(29) 



At a late time, the energy density of the curvaton and radiation from the inflaton decay at 
H ~ r CT are given by 



l ^ 4 f Civ 



and 



7T" 



Pr = -zz9*T R — - 



OR 



30^ - V«r a 

The energy density ratio of o to radiation at late time H ~ r CT is evaluated as 

.4 \ 1/4 



" 36MJ, V90M|r2 
For the case of v < cr os < cr v , we obtain 



YIMl (3(j2 s - u 2 ) VSOM^r 2 



(30) 



(31) 



(32) 



(33) 



IV. POWER SPECTRUM AND NON-GAUSSIANITY 



The curvaton is light during inflation with Eq. (Tl9l and thus has a Gaussian quantum 
fluctuation with the amplitude 5a* ~ H*/(2ix). The curvaton field value at the onset of 
its quadratic oscillation is some function of that at the onset of quartic oscillation o"2 OS = 
<T2os(c"os)- In addition we assume the field value at t os is same as that at horizon exit t*, 



Thus we can be expand cr 2o s around the homogeneous part a 2o s, 
where the prime denotes the derivative with respect to 0*. 



(34) 



(35) 



The curvature pertur 



using 8N formalism 29 



Dation due to the curvaton density fluctuation is easily calculated 



30[. The nonlinear curvature perturbation of the curvaton field on 



the uniform curvaton density hypersurface is given by 30|, |52|, [53] 



( a = 5N+- 



p(t,x) 



/■'II 



dp 

(t) P + P' 



(36) 
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where 5N is the perturbed expansion, p and p are the local density and pressure of the 
curvaton respectively. For the oscillating curvaton field in the expanding Universe, the 
pressure and the energy density is related by p = wp with w = when quadratic term 
dominates and w — 1/3 when quartic term dominates. In our case, the curvaton energy 
density at late times (t ^> t v ) is given by 

and it can be expanded around the background value a2 OS (t) using Eq. (135|) . 
From this, we can find the curvature perturbation of the curvaton field as 

Cx = Crl + \Co2 + ^Ca3 + ■ ■ ■ , (38) 



where each terms are expressed as 34] 



Cri = P^S**, (39) 

<JCT 2os 

U = -\ ( i - C 2 i - (40) 

U = 1(2- 3^^ + ) C\ - (41) 

4 V <T 2os <T 2os / 

On this uniform density surface at the curvaton decay time to or H = T^, we have 

Pr(t D )e 4 ^-0 + p CT (t D )e 3 ^-0 = Ptot (t D ), (42) 

with the radiation perturbation C- originated from the inflaton 0. Then, the curvature 
perturbation after the curvaton decay can be expressed as 

c = Ci + k 2 + k 3 + ..., ( 43 ) 

2 o 



where 



Cl=(l-i*)Crl+i*Crl, 

C 2 =(l - R)(r2 + RC.2 + R(i -R)(3 + R) (Cri - Ci) 2 , 

C 3 =(1 - R)(r3 + R(*3 + 3i2(l - R)(3 + R) (Crl ~ Crl) (Cr2 - Cr 2 ) 

+R(l -R)(3 + R)(-3 + 4R + 3R 2 ) (Cri - Cxi) 3 , 



(44) 



and 



Apr + 3p a 
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at t = tr>. Furthermore it is natural to assume ( r <C ( a and ( r is almost Gaussian so that 
only ( r i is non-negligible in the expansion of ( r . Thus in the above we can approximate 

Crl - Cxi ~Cal and Cr2 ~ Cx2 ^ ~Ca2- 

The power spectrum is obtained as 

V c = {l-RfV Cr + R 2 V^. (46) 

by using Eqs. (I3"8j) and f)43|) . In Eq. ()46p . the spectrum of radiation and the curvaton is 
given by 



V Hl_ 

(47) 



^ 2 ^2< s 



4vr 2 V 3a 2l 



with 



2 V V J ' H 2 

We defined r as the ratio of the contribution to the linear perturbation of ( from the curvaton 
to that from the inflaton, i.e. 

2-n D 2 /r. ' \ 2 



R 2 V Cr R 



:2e 



f^os) M 2 p _ (49) 



(i-Ryr (r (i-i?) 2 \3a 2o J p - 

In the limit of f — >■ oo, Eq. (T4"6"j) becomes that of simple curvaton scenario neglecting the 
inflaton contribution and in the opposite limit f — > the power spectrum has dominant con- 
tribution from that of the inflaton. Although for both cases the observable non-Gaussianity 
is possible [l^], the constraint from tensor-to-scalar ratio disfavors the small f region as we 
will see later. 

The scalar spectral index is given by 

n s = l-2e+^, (50) 



and we find 2e ~ 0.04 from the WMAP data, n s ~ 0.96. With this value of e, Eq. flU 
and Eq. ( 14T1) relate if* with the primordial power spectrum for the given initial value of <t os . 
Therefore we can obtain the each contribution from the curvaton and the inflaton separately 
and thus f in our scenario. 
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The tensor perturbation (gravitational wave) is also generated during inflation 24j . The 
tensor to scalar ratio rx is given by 



16e 



r T 



(51) 



P c (l-R) 2 (l+r) 

If the amplitude is large enough, the gravitational wave background is detectable through 
the measurement of the B-mode polarization in the cosmic microwave background (CMB) 
anisotropy by Planck j^J as well as the direct detection by future interferometers such as 



DECIGO 



26) 



We obtain the nonlinearity parameters 

5 f 2 



NL 



6 (1 + f) 2 



3 + A 2 
R 



2-R 



25 f 



,~3 



9nl 



54(1 + f) 3 



9 + 9A 2 + A 3 18 + 6A 2 



- 4 - 3A 2 + 10R + 3R 2 



(52) 



(53) 



R 2 R 

Here we have assumed that ( r is Gaussian so that only £ rl is non-zero and Cri Co-i which 
is true in the curvaton scenario when a* <C Mp. 

The WMAP data on the power spectrum of scalar and tensor perturbation 



V c ~ 2.4 x 10~ 9 , r T < 0.36, 



as well as the local type non-linearity parameter 



-10 < / NL < 74, 



(54) 



(55) 



constrain possible values of if*, cr* and R. 



A. small initial expectation value 

For a small initial amplitude of the curvaton field v < < cr v , the oscillation starts when 
the quadratic term dominates. In this case, from Eq. (1T8"]) with a os = a*, we find that 



a. 



2os 



a. 



(Tn 



(m a /2V\) 3/4 ( 


' 9<7 os (cr os - u)\ 


(3a 2 s - v 2 )^ \ 


, 4 (3a 2 os - v 2 ) ) 
( 9 9. 

/ 27 891a 2 s + (a os - 


{3a 2 s _ v 2yi/s 
(m CT /2 v / A) 3 / 4 


(3^-^)11/8 


V 4 + 16 (3a 2 s - 



_// , 297 (a os - v)a os 

2os 



r 3 



16 (3a 2 os -v 2 ) 
v)a os 16929 (cx os - v)a d 



(56) 



! I - y2) g4 _ ^2 
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The curvature perturbation of the curvaton ( a and the corresponding non-linearity parame- 
ters are evaluated from this. Note that there are additional factors coming from the depen- 
dence on H 2os compared to the simple curvaton model with quadratic potential. One thing 
to note is that A 2 , defined in Eq. f ]4"0"]) . becomes negative when a os > 1.8v, which changes 
the sign of /nl in Eq. (152 p for small R. The non-trivial behavior from this is shown in the 
lower part of Figs. H] and |5j 

For the pure quadratic potential limit <r s — v, the expression of ( a is reduced to 

c=-— -l/— y+-f— y. (57) 



3 cr* — V 3 \ cr* — v J 9 \ a* 
The density perturbation of radiation after the curvaton decay is 



(58) 



The nonlinearity parameters are given for v < a* < a v , from Eqs. ( )52j) . (1531 and (1381 . 



r \ 2 5 / 3 „ \ 3 



f "^{— f ) l\2R- 2 - R ) > -2 

( f \ 3 25 / 9 1 9 
9 -=(tT?) 54(-S + 2 + 10fl + M 

The sizable large /nl ~ 100 is obtained with a small ratio R ~ 10" 



(59) 



B. large initial expectation value 

1. evolution of perturbations 

Next, we consider a large initial amplitude of the curvaton field a* ^> v. At the early 

stage of oscillation, the field evolution is due to the quartic potential and highly nonlinear. 

For our model, the corresponding quantity <t 2o s is analytically related to a* using Eqs. f lT5|) 

and ( 154"]) . ignoring the B part. Then, we obtain 

1 da 2os 3 
c"2os da* 4a* 

(a') 2 ~ 3' [W) 



'2os 

rr 2 a'" 5 

°2os°2os 



^2os) 3 9' 
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for a high reheating temperature case, and 

1 da 2os 1 
o"2os dcr* 2a* 



ff2osff2os ~ -1 (61) 
3, 



2os/ 



for a low reheating temperature case. The prime denotes the derivative with respect to a*. 
In addition, concerning with the part, because of the high nonlinearity, the fluctuation 
5a also undergoes nontrivial evolution. The equation of motion for 5a of superhorizon scale 
(k <C Ha) is given by 

5a + 3H5a + X(3a 2 - v 2 )5a = 0. (62) 

During a ~ 0, 5a has effectively the negative mass. This tachyonic instability leads to 
significant amplification of the fluctuation 5a, in some cases that the initial value <j os cor- 



responds to the transition of the VEV from — v to v in Fig. [TJ |33|] • Figure [3] shows the 
amplification and evolution of the field fluctuation 5a and the density ( a for some a os s. For 
cases in which the field a stays near the origin longer, the amplification is sizable and, with 
Eq. ( JIB"]) , roughly estimated as 

rp _ 5a2os ^ e VxvAt < e Vxvt v „ e ^ (63) 

OCT* 

for a high reheating case, as seen in the middle row in Fig. |3]for <r os = 25v. Here At denotes 
the period during the tachyonic instability works. The final fluctuation after quadratic 
oscillation starts is, with the amplification factor T, given by 

H 2os = T^fc, (64) 

However, this enhancement occurs only for limited conditions of a os near the VEV transition 
intial expectation value. Thus, from now on, we consider cases without this enhancement 
and these enhanced modes will be studied in future works. Then, we obtain the curvature 
perturbation of the curvaton field 



CrU t „ = \— ~ T ( — ] + 7 I — ) , for hi g h T Ri 



2 (7* 



5a* \ 2 




( 5a* 


o-* J 




\ o* 


5a* \ 2 


I 


( 5a* 


(J* ) 




\ o* 



^4^-6^ + 9^J ' f0rl0WTR ' 
This is conserved until the curvaton decay for t^$>t v , as seen in figure. [31 
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(65) 
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FIG. 3: The evolution of a (left), 5a (center), and £ CT (right). The upper (middle, lower) row cor- 
responds to the results for <7 os = 5 (25, 50) v. Notice that the vertical axis of the right-middle figure 
is logarithmic scale. Here, we assume 5a* = H*/(2ir) with H* = 10~ 3 v. The green (purple) line in 
the right figures expresses the analytic formula of £ CT at t os (t v ) without including the amplification 
effect. These show that the error of the analytic formula is just about 0(10)% unless the tachyonic 
instability is induced as for <7 os = 5u, 50t>. 

2. observables 



Since the total curvature perturbation of radiation after the curvaton decay is conserved, 
it can be calculated at the time of the curvaton decay. For the case of double well poten- 
tial, i.e. initially quartic term dominates and the mass term becomes dominant before the 
curvaton decays, the primordial curvature perturbation is obtained, from Eqs. and ( 165]) 
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with a os = a*, as 

C = (i 



A \ a 



+ 



48 









:) 
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+ 2 


^2 


~ R 





5 a* 

a* 

10R + 3R 2 ) R 



and 



Jiff, 



5 a* 
a* 



1 (-2 

18 v 7 W* 



(66) 



for high T R , 



1 / A \ ^ 

+ 162 ( 5 + 10jR + 3jR ^ ^ W ) ' f ° r l0W Tr ' 



(67) 



where we assume that the perturbation of radiation, which has the origin from the inflaton 
field ( r , is Gaussian and these higher order contributions are negligible, compared to those 
of the curvaton. Here R is evaluated when the curvaton decay as 



R 



3p c 



4p r + 3p a 



3r 



H=T„ 



4 + 3r 



(68) 



H=T„ 



where p r and p a are the energy densities of radiation and the curvaton respectively and r is 
the ratio of them, r = p a j ' p r . Note that ( r can be comparable to R( a with small R, which 
becomes the general mixed inflaton-curvaton scenario 32]. 

The nonlinearity parameters are given for a large initial amplitude, a* ^> a v , for high Tr 
from Eq. (|66|) as 

2 



NL 



9nl 



r 



1 + rJ 6 \R 

3 



'l(L- 2 -R)>-2 



r 



1 + r 

and for a low Tr from Eq. fl6T|) as 



54 1 R 2 R 



3R 



(69) 



for high Tr, 



/] 



NL 



fl'NL 



r 



1 + f/ 6 

3 



/?) 



r 



r 



25 

— (5 + 10/? + 3/? 2 ) 
54 v ' 



(70) 



for low Tr, 



The sizable large /nl ~ 100 is obtained with a small ratio R ~ 10~ 2 for a high Tr. However it 
is impossible to have such a large non-Gaussianity for a low Tr since there are cancellations 
in the coefficients of inverse of -R-terms. In this region, the curvaton starts the quartic 
oscillation, when the Universe is dominated by the oscillating inflaton field. 
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FIG. 4: Left: The contour plots of R and r T . Here we used v = 10 15 GeV and T R = 10 12 GeV. 
Right: The contour plot of /nl in the same plane of the left. For the sign of /nl refer to figured 
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FIG. 5: Left: The contour plots of R and ry. Here we used v = 10 15 GeV and Tr = 10 14 GeV. 
Right: The contour plot of /nl m the same plane of the left. For the sign of /nl refer to figure 

In Figs H] and [5j the contour plots of gravitational wave background and R (left window) 
/nl (right window) are shown with observational constraints for the decay rate given by 
Eq. ([6]) in the plane parameters of A and cr*. Here we have fixed VEVs v = 10 15 GeV, 
and the reheating temperature T R = 10 12 GeV (figure HJ) and T R = 10 14 GeV (figure [5]) 
separately. The observational constraints include the tensor-to-scalar ratio = 0.1 and 
10~ 3 for the expected sensitivity of B-mode detection by Planck and future instruments 
such as CMBPol [54|, respectively. For given <7 os and A, imposing Eq. ( )54|) on Eq. ( 1461) 
determines if* with e ~ 0.02 from Eq. (150j) . 

As can be seen in the figures, the large a* or the large A region is excluded by the null 
detection of gravitational wave. The region with f < 1 belongs to this excluded region, 
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FIG. 6: The plots for / NL (left) with the fixed A. Here v = 10 15 GeV and T R = 10 12 GeV (left) or 
T R = 10 14 GeV (right). 

which means that in the allowed region the power spectrum is dominated by that from the 
curvaton in our scenario. Above the diagonal line, written above (curvaton dominated), 
the curvaton energy dominates the Universe before it decays. In the limit of curvaton 
domination (R = 1) the non-linearity parameter becomes /nl = —5/3 or —5/2 for high Tr 
and low Tr respectively. Below the line '(curvaton dominated)', R can be much smaller 
than 1 and thus there is a chance to obtain large non-Gaussianity for high Tr case. This 
region appears in figured] (Right window), as a wedge shape above <j os > a v . In this region, 
it is possible to generate large non-Gaussianity of the order of 20. On the other hand, for 
a larger initial amplitude the curvaton oscillation starts before reheating, corresponding to 
low- Tr, and the non-Gaussianity is significantly suppressed because of the cancellation as 
discussed after Eq. flTOj) . This happens in the '(low Tr)' region above the dashed diagonal 
line in figure |5j For Tr = 10 12 GeV (figure H]), all the drawn region corresponds to '(low Tr)' 
and the non-Gaussianity is small for <t os > a v . 

The interesting behavior happens for the initial amplitude v < er* < a v , where the 
oscillation of the curvaton starts in the quadratic term dominated potential. There are 
two regions depending on the sign of /nl, positive for v < cr* < 1.8v and negative for 
1.8v < a* < a v . This difference is due to the evolving effective mass in a double well 
potential, which is constant in the pure quadratic potential. In the figure [6] we show the 
plot of /nl depending on cr* for given A for Tr = 10 12 GeV (left) and 10 14 GeV (right) 
respectively. The positive /nl with the magnitude of the order of 100 is possible for a small 
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initial amplitude v < o* < 1.8v. 
V. SUMMARY 

We have studied the density perturbation generated by a curvaton whose potential is flat 
with small self-coupling in a double well type. We have used a large VEV, v > 10 15 GeV, to 
avoid the domain wall formation and a larger initial curvaton amplitude a* which is easily 
obtained in this flat potential. 

We have analyzed the cosmological evolution of the scalar field in this flat double well 
potential to see the viability of the field as a curvaton to generate the primordial density 
perturbation to explain the structure formation and the anisotropies in the CMB. With 
a large initial expectation value a os ^> v, the energy density decreases as initially oc a -4 
and at late time as oc a~ 3 . This is same as a massive curvaton with self- interaction with 
vanishing VEV. However, three crucial differences and features appear. One is the tachyonic 
amplification of the fluctuation by the negative mass squared at the origin, for a particular 
initial value a os which periodically appear in the parameter space. We have shown that, 
except for these tuned boundaries, the density perturbation and other nonlinear parameters 
are well approximated by analytic formula. The second is the suppression of non-Gaussianity 
even for very subdominant curvaton, if the reheating temperature after inflation is as low as 
not to satisfy Eq. ( I23p . The other is the non-trivial behavior of /nl even when the curvaton 
field starts oscillation trapped at one of the minima. We found a successful scenario needs 
the flat potential with small self-coupling of the field of the order of < 10 -10 for a reasonable 
reheating temperature. 

In conclusion, we have found the differences in massive self-interacting curvaton models 
with and without VEV. In addition, in a double well potential curvaton model, non-linear 
parameters can not be so large for a large initial field value and a low reheating temperature. 
Therefore, if both large nonlinearity and B-mode polarization will be detected, the potential 
of the curvaton or thermal history of the early Universe will be constrained. 
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